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XSSQRT 3; Floating Point Square Root ele = 6-$ 4 AX/VMS ™ v04-00 P 
ott ° ” need aces P=] ope 04: 98: ¥ MTHRTL. eae eooGRy MAR; 1 — 


~TITLE UVXSSQRT 5 vegat gating Point Square Root routine 
. IDENT /1-016/ : rite:  ATHSORT.MAR EDIT JCW1016 


SARAH AAAAAAAAH AAA AEAAAAEAAAAAAAAAAAAAEARAAARAAAAAAAAAARAAAAAAEAAAAARAAAAEAAEEEREEE 


COPYRIGHT (c) 1978, 1980, 1982, 1984 B 
DIGITAL EQUIPMENT CORPORATION, MAYNARD. MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


® 
* 
® 
® 
® 
THIS SOFTWARE IS FURNISHED net A LICENSE AND MAY BE USED AND corer’ ® 
ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE * 
INCLUSION OF THE ABOVE 3 gill NOTICE. THIS SOFTWARE OR ANY OTHER * 
COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY ° 
TRANSFERRED. e 
® 
Hs 3 INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE * 
AND SHOULD NOT * 
® 
® 
® 
® 
® 
® 
* 


BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 
CORPORATION. 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


COOoooooooooooooooooooo 


SOCOCCCOOCOOOOOOOOOOOOCoOO 
OOOO COoOOoCooCooooooo 
DOOCOCCCOCOCOCCOOCOCCoO: 


. . 7 . .* . ’ 
eaeeeeneeeeaneeeaeeeenneaeanne 


IRR RRAASAAAASAAAAALASAAAAAS ALARA LARASESASARRR RARER ARERR RASA SA SERRA RAR AAAS SD SO 


0000 
FACILITY: MATH LIBRARY 
"ABSTRACT: 


MTHS$SQRT is a function which returns the floating point square root 
of its single precision floating point argument. The call is standard 
ATHSSORT R erence. 

MTHS$SQRT_R3 is a special routine which is the same as MTH$SQRT except 
a faster non-standard JSB call is used with the argument in RO and no 
registers are saved. 


VERSION: 01 


HISTORY: 
AUTHOR: 
Peter Yuo, 15-Oct-76: Version 01 


MODIFIED BY: 
1-3 Peter Yuo, sfonet7 


1-2 Peter Yuo, 31-May-7 
1-016 Jeffrey ¢. Wiener bon eb-83 


OOOO QOCoCoooooooooo 


SOSCSCSCOCSOSOOSOOOSOSOSOOOOSOOSOoOoSoSo 


CSOSCSCCSOOCOSOSCOOOCOSOOOOoOOOOoOooo 
WRI O ODNA UNE WIN 9 OO NA UE WIN | O OO NIA UE WIN @ O OO NAU EW OOONAU EWN 


i paieaeadie Reet vist Vung ee apie on Sete Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Se Ge Ge Ge Ge Ge Ge 


COOOSOOCSOOOSOSOSOOSOOSOOOOOOOSSSoSoo 
PVP BP BB BBE BE NANNING RININIPUPONPINININ 2 PO OO ee 


oooooooo 


So 
So 


-C 


] 
Hi 
XSSQRT ; Floating Point Square Root Ra SEP-1984 02:08:14 VAX/VMS Macro V04-00 Page 
ote HISTORY Y : Detailed Current Edit History ¥? p= 198 95:98:38 MTHRTL.SRC REI VXSQRT.MAR; 1 . 3) | 
.SBITL HISTORY ; Detailed Current Edit History | 
| 
8 ; ALGORITHMIC DIFFERENCES FROM FP-11/C ROUTINE: none | 
¢ ; Edit History for Version 01 of MTHSSQRT | 
66 3 Se Code saving after code review 
0 65 ; 01-2 ROTL shift in on afte into highest bit. Use ASHL instead. 
+4 64 ; ADOL instcycs ton after ADJUST has sabes: 6 changed into ADDW to prevent 
0 65 ; over flow if R1<31:16> = FFFF on 
. 66 ; §1-3 Finish error ened tae 10-June-1 
67 ; O1- he aga chan .* te MTH SToNaL. 
00 68 ; MTH$_... change Mo» 
000 +4 ; Changed error Sendl lan mechanism. Put error result in RO before 
B88 0; calling MTHSSSIGNAL in order to allow user modify error result. 
00 71 3 p1-6 Return -0.0 on negative arg. TNH 20-Dec-77 
44 i : pi Edit in Rich Lary’s code bums. JSB rousing f now, _R3. JMT 19-Jan-78 
44 73; 01-9 Move .ENTRY symbol to module header. TNH 1 aug 
0000 74 ; 1-010 = Put version number 2 stentere format: # at Wes ed 
0000 3 3 numbers. Also updat Mae lane notice. gs 6 -Nov-78 
0000 76; 1-011 - Change, nth SQUROONEG to tatu Q UROONEG. 18s 0 7=DEC-78 
0000 77 3 1- \¢ ~ o the PSECT directive. WBS ¢ 
0000 78 ; 1-013 - Declare aoenreaha. SBL 17-May-19 
0000 79 ; 1-014 = Move MTHS$SQRT_R2 to Sepgrate sogule (MTHSQRTR2.MAR) and 
0000 80 ; replace with ATHSSQRT JAW go-8 ope 1979. 
+44 81 ; 1-015 = Changed W* to G* in call to MTHS$SIGNAL RNH 09-Sept-1981 
000 Hf 3; 1-016 = Changed all ror to D_floating to G_floating. This was 
0000 83 ; done because & floating, not D floating, Ts to be used as 
0000 84 ; back-up to F_floating. JCW 22-FEB-8 


—co 
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Point Square Root ro sles 16-SEP-1984 :08:14 VAX/VMS Macro v04-00 P 
Bus’: Bes turatire Bare st Modul Gasepatobe $fi58i56 PATANTE eaeSuvyccns mans1 PM Cd, 


-SBTTL DECLARATIONS ; Declarative Part of Module 


XSSQRT ; Floa 
oe BECLAR 


: 

: ; INCLUDE FILES: 
91° 

38 


: ; EXTERNAL SYMBOLS: 
95 ° .DSABL GB 


L 
~EXTRN MTHSK_SQUROONEG 
eEXTRN MTHSSSIGNAL 


: EQUATED SYMBOLS: 


0000400C ACMASK = “M<IV, R2, R3> 3; register save mask and IV enable 
3; MACROS: none 
> PSECT DECLARATIONS: 
0000 -PSECT _MTHSCODE PIC,SHR,LONG, EXE ,NOWRT 


3 program section for math routines 
OWN STORAGE: none 
CONSTANTS: 


Ok kd dt 
RMNMNN 2 2 2 OO SS S“QOOOOOOOCOCOOOOOOOO 
FWN OOO NAU EWN 0 DOONAN EWN 0 OONO 


0 
0 
0 

0000 

0000 : 

44 3; Constants A and B chosen for k = odd 
13CD5FD4 $009 ; LF_ODD_A_£63_ = “X13CD5FD4 
3C€4A2018 44 : LF_ODD_B_EM63 = *x3C4A2018 

444 : Constants A and B chosen for k = even 
F61A4015 $000 : LF_EVEN_A : “XF61A4015 
48231FD7 0000 LF_EVEN_B_EM64 = *X4B231FD7 


-cC 


SOOoooooooo 


OUTPUT PARAMETERS: 


VALUE: floating square root of the argument 


IMPLICIT OUTPUTS: none 


5 
XSSQRT ; Floating Point Square Root rou « 16-SEP-1984 02:08:14 VAX/VMS Macro v04-00 Pa 4 
=016 WTHSSORT 9~ Standard Single Provision FL Gasepe1obe $ii58,26 EATARTE eaedovKcoRsomarn;1 Pe! chy 
! § -SBTTL MTHSSQRT = Standard Single Precision Floating SQRT 
158 s+ 
: 9 ; FUNCTIONAL DESCRIPTION: 
! § : - single precision floating point function 
99 , 3 SQRT(X) is computed using the following approximation technique: 
0 0 : § : If X <= 0, error. Let X = iXi. 
p 9 : 8 ; Let X = 2eeK * F where F is the fractional part. 

909 140 : If K = even, X = 2e0(2P) * F 
99 141 : SQRT(X) = 2eeP * SORT(F), 1/2 =< F< 1 
000 128 : If K = odd, X = 2e2(2P+1) * F = 2e2(2P42) * (F/2) 

900 146 : SQRT(X) = 2ee(P+1) ® SQRT(F/2), 1/4 =< F/2 < 1/2. 

9000 146 : Let F' = AeF +B, 

000 147; A = 0.453730314(octal), 

4444 148 ; B = 0.327226214(octal), for K = even. 
000 149 ; = Ae(F/2) + B, 

000 150; A= 8 Seer eeioctal). 

000 1517; B = 0.230170444(octal), for K = odd. 
000 13 ; and 

0000 155; K' =P, for K = even 

0000 154; =P +1 for K = odd. 

0000 155; 

0000 156; Let YO = 2**k' * F' as a staight Line angres tastes wthin the 
000 157; given interval using coefficients A and B which minimize the 
44 138 3 absolute error at the midpoint and endpoint. 

0000 169 : Starting with YO, two Newton-Raphson iterations are performed. 
00 168 : y{n¢1) = (1/2) * ( YEn] + X/YLEn)) 

00 163 ; 
000 164 ; The relative error is < 10**-8, 
000 165; 
000 166 ; CALLING SEQUENCE: 
4 167 ; 
00 168 3 sart.wf.v = MTHSSQRT(x.rf.r) 
33 1%6 : INPUT PARAMETERS: 
00000004 00 76 LONG = 4 3; define Longuord multiplier 
00000004 We x = 1 * LONG ; Contents of x is the argument 
175 ; IMPLICIT INPUTS: none 
176 
1) 
135 
18) 
181 
182 


yurssant Arf bee ng_ Point Square Root rou tine” tT hs 96:08:16 AX/VMS Macro V04-00 Page 5 


ndard Single preci sion + =SEP-19 4 MTHRTL.S REI VXSQRT.MAR; 1 (4) 
COMPLETION CODES: none 
SIDE EFFECTS: 


Signals: MTHS_SQUROONEG if X < a -0 Ay reserved ooerand | in RO (copied to 
gnal mechanism verter ake’ LM gh a LIBSSIGI 

Associated message is: ROOT OF NEGATIVE VALUE"’. Result is reserved 

operand -0.9 unless a it. Socekiaa (or any) error handler changes CHF SL -MCH_RO/R1 


NOTE: This procedure disables floating point underflow, enables integer 
overflow, causes no floating overflow or other arithmetic traps, and 
preserves enables across the call. 


400C eENTRY MTHSSQRT, ACMASK 3; standard call-by-reference entry 

3; disable DV (and FU), enable IV 
MTHSFLAG_JACKET ; flag that this is a jacket procedure in 
6D O00000000°GF 9E MOVAB G“MTHSSJACKET_HND, (FP) 

3 pat anne? address to jacket 

3; handler 


8 3 case of an error in special routine 
0 MOVF ax(AP), RO ; RO = arg 

04 BSBB MTHSSQRT _R3 3; call specail SQRT rountine 

05 RET 3 return - result in RO 


at 


eH 


ng Point Square Root feattieg 16-SEP-1984 AX/VMS Macro V04-00 Pa 
as °- jsB Sar! routine -SEP-1984 t 9% 98: ¥ MTHRTL.SRCJUVXSQRT.MAR; 1 - (8) 


yerssont 

-016 
-SBTTL MTHSSQRT_R3 = JSB SORT routine 

JSB SQRT = used by the standard, and directly. 


CALLING SEQUENCE: 
save anything in RO:R2 
MOVE 


JSB MTH$saRT_R3 
return with result in RO 


; input in RO 


Note: This routine is written to pe ony integer overflows, floating overflows, 
floating underflows or divide by 0 conditions, whether enabled or not. 


saath * USED: 

- Floating argument then result 

RY - X saved for use during iteration 
R2 = scratch 


sme’ 


JSB routine for SQRT 


PEO MOO MS FW DUNNO OOOOOOCOOOOCOOOOCOOCOCSCO I 


" 
3 
3 
4 
4 
4 
4 
4 
4 
& 
4 
4 
4 


1 
1 
18 
14 
12 
i 
18 
+ 
4 
7% 
51 50 50 $ RO, R1 ; test sign of X and save it in R1. 
62. «+15 $ BLEQ ZERO_NEG ; branch to ZERO_NEG if X =< 0 
0:x>0 
! bos: 
S20 : ; oe ee RO, R2 : isolate low 16 bits +4 oup.> tract) in R 
5 94 4 CLRB p 3: R2 now has i and aoe ? que b 
50 5 AA 5 BICW R2, RO ; clear sign an hett 7 exp bits P 
5 95 6 TSTB RO ; check low bit of exp 2 
10 «618 7 BGEQ EVEN 3 and branch I 
50 13CD5FD4 B&F 44 8 MULF #LF_ODD_A_E63, RO 3; add 64 (half of bias) to (exponent-2) Cc 
.) 3 and start approximation calc P 
50 3C4A2018 8F 40 0 ADDF #LF_ ODD_B_EM63, RO ; RO = (first approx) * 2%*-64 s 
13° «11 1 BRB ADJOST ; go adjust 4 
: EVEN: P 
8F Ad 4 ADOW *x2000, no : exp is 0 - make it 64 (2**-64) for legalit Cc 
8F 44 00 5 MULF At EVEN RO A 
8F 40 00 § ADDF #LF-EVEN “8° "EM64, RO ; RO = (first approx) * 2**-64 
ADJUST: T 
9C 8 ROTL #31, R2, R2 3 divide R2 Cexpsbias) by 2, 3 
0 9 : giving (exp/2+64) T 
AO 004 50 ADDW R2, RO 3 inser exp/e in first approx and 3 
re 2] 3 re-bias it 1 
rt 8 first iteration - single precision is sufficient 
a ne a tae 
A2 05¢ 38 SUBW BiaG. RO RO = 11 = (1/2)(¥0 + x/Y0) 
be 8 : no overflow possible 
057 60 : second iteration, do in double precision to get truncated( rather than 
037 g3 : rounded) result. 
57 e§ | oe3 CLRL . ; Lower part (xX) = 0 
057 64 33; DIvD RO, R1 3; divide Y1 into X with low-order 


a _— 


Mm 5 
Point Square Root routine 16-SEP-1984 02:08:14 VAX/VMS Macro v04-00 Page 7 
As - jsB SORT routine pane} 94:98:38 MTHRTL.SRCJUVXSQRT.MAR; 1 . (5) 
; 32 bits of Y1 garbage. This doesn't 
; effect accuracy, since Y1 inacurate 
3 anyway. 


ivide Y1 into X with low-order 53-24= 
9 bits of Y1 garbage. This doesn't 
effect accuracy, since Y1 inacurate 


$36 Using oy instead we have 
338 DIVG RO ,R1 


RO goes to signal mechanism vector 

(CHF SL_MCH_RO/7R1) so error handler 

can modify the result. 

signal error and use real user's PC 
independent of CALL vs JS 

return - RO restored from CHFSL_MCH_RO/R1 


00000000°GF 02 FB 


CALLS #2, G*MTHS$SIGNAL 
RSB 
- END 


ti 
RT 
$3 
: 
49 
i anyway. 
52 1 99FD 74 CVTFG RI, R2 3 convert and copy X into R2/R3 
1 04 75 CLRL R1 ; clear low part of Y1 
43 0 99FD 5 6 eyes RO, -(SP) 
5 E 46FD 06 f f DIvG2 (SP)+, R2 ; divide ¥1 into x 
06 i : Need to CVTIGF without rounding. To do this clear the ey rc bit. 
3 80 ; The rounding bit is bit 12 of the 3rd word(word 2 starting from 0). 
53 FFF BF AA Boe Ht . BICW2 #*XIFFF, R3 3; Clear the ane bit so that 
52 52 33FD One HH k CVTGF R2, R2 ; CVIGF truncates instead of rounds. 
B98 HH : Back to original code. 
50 52 40 $08 a : ADDF _—R2,_RO : RO = 1 + higher_part(x/Y1) 
50 0080 8F Ag 007 88 SUBW #*x80, RO 3; RO = SQRT(X) = (7/2) (Y1 + X/Y1) 
05 007 89 SQRTX: RSB : return, RO = result 
007 90 
007 91; x =< 0 
007 3 . 
007 93 2ERO_NEG: 
FD 13 007 94 BEQL SQRTX return with RO = result = 0 
6— dD 007 95 PUSHL (SP) return PC from JSB routine 
7E OO'8SF 9A th 96 MOVZBL #MTHS$K_SQUROONEG, -(SP) condition value 
50 01 OF 78 3 ze ASHL #15, #T, RO RO = result = reserved operand -0.0 
6 8 99 
008 00 
008 4 
b08a 308 
0 0 
0 0 


Be ere 8 Osos sO rimismiim Se nen OW aoe a od ee 3s 
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; Floating Point Square Root routine 


UVXSSQRT 
Symbol table 


ACMASK = 4006 
ADJUST 44 R 01 
EVE 31k 1 
LF_EVEN_A = FEIAKOI9 
LF EVENT BE M64 4p231FD 
LF“ODD_K_€63 CDSFD4 
LF“ODD"B-EM63 = 3C4A2018 

G 0000004 


= 
MTHSSJACKET_HND tteeeeee 
HS$SSIGNAL eeeeeree 


MTHSK SQUROONEG eeeeeeee Xx 0 
MTHSSORT 9000000 RG 0 
MTHSSQRT_R3 00010 RG 

0 R 8 


SQRTX 00000076 R 


=< 


The working get Limit was 900 pages. 

4312 byies (9 pages) of virtual memory were used to buffer the intermediate code. 

There were 10 pages of symbol table space allocated to hold 17 non-local and 0 local symbols. 
5 source Lines were read in Pass 1, producing 11 object records in Pass 2. 

1 page of virtual memory was used to define 1 macro. 


$eoeceeeeeeoeoccosocoscoscos 


: Macro Library statistics : 


Macro Library name Macros defined 


_$255$DUA28: (SYSLIBISTARLET.MLB; 2 


"SSEP=1986 99:29:26 EATHNTL.SRESUVNSGRY MAR; 1 


¥ = 00000004 
ZERO_NEG 00000077 R 01 
pees er een eer cawaece + 
! Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
- «ABS. 00000000 ( 0.) 00 ¢ OO.) NOPIC USR CON ABS LCL NOSHR NO 
_MTHSCODE 00000088 ¢ 139.) O1¢ 1.) PIC USR CON REL LCL SHR 
ee + 
! Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 29 $88 08-1) 00:00:00.68 
Command processing 107 Bo: 00.64 Se 
Pass 82 0:00:00.81 00:00:03.90 
Symbol table sort - 1 -f 0:00:00.00 
Pass 2 6 :00: 8 0:00: fet] 
Symbol table output 0:00:00. 0:00: 6) 
Psect synopsis output 3: :0 8 8: 0:00. ‘ 
Cross-reference output :00:00. 0:00:00.0 
Assembler run totals 29 00:00:02.2 00:00:09.93 


Ee — = -_ 


v04-00 P 
) age (3 


EXE NORD NOWRT NOVEC BYTE 


} 
RD NOWRT NOVEC LONG 
| 
| 
| 
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UVXSSRT :; Floating Point Square Root renine 
VAX-11 Macro Run Statistics 


0 GETS were required to define 0 macros. 
There were no errors, warnings or information messages. 


"SSEb19RE 99:35:28 


AX/V 
MTHR 


TL.SRED 


em eS NS a 


o v04-00 
UVXSQRT.MAR;1 


re, 


MACRO/ENABLE=SUPPRESSION/D1 SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:UVXSQRT/OBJ=OBJ$:UVXSQRT MSRC$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) +MSRCS: 


NT CORPORATION 
PROPRIETARY 


AH-BT13A-SE 
VAX/VMS V4.0 


| 2b 


